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Abstract 

After a short introduction on particle candidates for dark matter within possible 
extensions of the standard model, we concentrate on Weakly Interacting Massive 
Particles, and on one of their most interesting physical realizations: the neutralino. 
We analyze how detectability of relic neutralinos by direct and indirect means is 
related to their local and cosmological densities; we use simple general arguments 
to discusss different scenarios where relic neutralinos make up the dominant bulk 
of dark matter or only a small fraction of it. Our general arguments are further 
corroborated by specific numerical results. We show to which extent the present 
experiments of direct searches for WIMPs, when interpreted in terms of relic neu- 
tralinos, probe interesting regions of the supersymmetric parameter space. Our 
analysis is performed in a number of different supersymmetric schemes. 



1 Introduction 



Evidence for existence of dark matter and dark energy in the Universe is provided 
by a host of observational data: properties of galactic halos and clusters of galaxies, 
large scale structures, cosmic microwave background, high-redshift supernovae SNe 
la. As far as dark matter is concerned, a favorite range for Vtmh^ {^m being the 
matter density divided by the critical density and h the present-day value of the 
Hubble constant in units of km ■ s~^ ■ Mpc~^) may be set as: 0.05 < flmh"^ < 
0.3. Notice that the most recent determinations of cosmo logical parameters 
appear to pin down the matter relic abundance to a narrower range 0.08 < Vtmh? < 
0.21; however, some caution in taking this range too rigidly is advisable, since some 
determinations of cosmological parameters are still subject to fluctuations. We point 
out that in the present paper we do not restrict ourselves to any particular interval 
of VLmh?] only some features of Figs. |^|, depend on the actual value employed for 
the minimum amount of matter necessary to reproduce the halo properties. 

2 Particle candidates for dark matter 

Various possible extensions of the standard model, envisaged in particle physics for 
reasons quite independent of any cosmological motivation, offer, as an extra bonus, 
a great variety of particle candidates for dark matter. 

A very natural candidate is the neutrino: both a light neutrino, with a 
mass iriy < 1 eV, or a heavy one, with a mass iriu > 100 GeV, might be viable 
and interesting candidates. The light neutrino will be briefly discussed in the next 
section. 

What appears to be the most natural (though certainly not the unique) 
solution of the strong CP-violation, the axion 0), would also be an appealing 
candidate for dark matter; in fact, though already constrained by a number of 
observational data, the axion still offers interesting cosmological perspectives. Lack 
of space does not allow here a presentation of the current situation; we refer to 
Ref. 0) for an overview of this subject. 

One of the most favorite candidates for dark matter is the neutralino, 
which in many supersymmetric schemes turns out to be the Lightest Supersym- 
metric Particle (LPS). Supersymmetry is known to be motivated by a number of 
strong theoretical arguments, though it is not yet supported by experimental evi- 
dence, except for a few possible hints: unification of coupling constants at the grand 
unification scale Q/, Higgs LEP events BSE) , muon anomalous magnetic mo- 



0) (the two latter points will be briefly discussed in Sect. 5). If 
supersymmetry exists in Nature and the R-parity is conserved, the LSP would be 
stable. Depending on the susy-breaking mechanism and on the sectors of the susy 
parameter space, the LSP may be one of many options: sleptons, squarks, gluino, 
gravitino, neutralino, axino (if the axion exists). When weakly interacting, the LSP 
has the prerequisites for being a good dark matter candidate. In the present note we 
discuss in detail the neutralino, whose detection rates may reach the level of current 
experimental sensitivities, both as a dominant or a subdominant component of dark 
matter. Other possible susy candidates are not discussed here. 
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3 Decoupling of particles from the primordial plasma 

In the standard scenario of the early Universe, relic particles originate from their 
decoupling from the primordial plasma, while in thermic equilibrium; the decoupling 
of any single species occurred when the interaction rate of that species with the 
thermic plasma became smaller than the expansion rate of the Universe (at the so- 
called freeze-out temperature). It is customary to deflne as hot (cold) candidates 
those particles that are relativistic (non relativistic), when they decoupled from the 
thermal medium. 

A typical case of hot candidates is provided by neutrinos with a small mass, 
rui, < 1 MeV. Their present-day rehc abundance turns out to be Quh"^ = I],^m,^/(93 
eV). Actually, we already have some indications of non-vanishing neutrino masses. 
The present experimental data may be formulated in terms of difference in squared 
masses, namely: Am^ ~ 10~^ — 10~^ eV^ from solar neutrinos (if MSW effect is 



at work) EH), Am^ ~ (2 — 5) x 10 eV^ from atmospheric neutrinos EIP, 0.2 
eV^ < Am^ < leV^ from the LSND experiment 113). 

Using the data on atmospheric neutrinos, a conservative estimate gives: 
^i/h^ ^ 5 X 10~^, i.e. a value of relic abundance which, although small, is already at 
the level of the average density of visible matter. Notice that Q^h"^ could be sizeably 
larger, in case of a degeneracy in neutrino masses. Unfortunately, at present, no 
experimental device is capable of measuring these low-energy neutrinos [T^ = 1.95 
K). Finally, we recall that light neutrinos cannot comprise the whole bulk of dark 
matter; indeed, should they constitute a dominant dark matter component, one 
would have difficulties in generating the present cosmological structures. 

A generic Weakly Interacting Massive Particle (WIMP) would behave as 
a cold dark matter candidate. Its relic abundance may be derived to be 



where (Tann and v are the WIMP pair-annihilation and the relative velocity, respec- 
tively. < aann V > dcnotcs the thermal average of the product {aann ■ v) integrated 
from the freeze-out temperature to the present-day one. What is remarkable is 
that for a number of particle candidates it is conceivable that Q^jMph"^, as given 
in Eq. (|l]), may be of order 0.1 — 1. Even more noticeable is the fact that for some 
candidates with fl^riMph'^ = O(O.l) the detection rates are at the level of current 
experimental sensitivities. This is the case of the neutralino, which will be discussed 
at some length in the following. However, it is obvious that detectability of any given 
candidate is of the utmost importance, also in the case of a relic particle which can 
only contribute for a fraction of the total dark matter density. We will examine this 
point in Sect. 4.3. 

One should be aware that also decoupling mechanisms different from the 
standard one are viable. An interesting example is provided by the low reheating- 
temperature scenario, discussed in Ref. 0). For instance, in this case the WIMP 
rehc abundance could be proportional to <Jann, at variance with Eq. (|l]). Thus, the 
relevant WIMP phenomenology could be quite different from the one depicted in 
the rest of this note. 

Dark matter particles may also be produced non-thermally (see, for in- 
stance, Ref. ED for a review and Ref. for a specific realization leading to relic 
neutralinos). 

4 Direct and indirect signals of WIMPs 

Different strategies for detecting the presence of WIMPs in our halo have been en- 
visaged. The present note is mostly devoted to an analysis of the discovery potential 
of the direct searches for WIMPs, though some considerations on the indirect means 
are also included. 

4.1 WIMP direct detection 

The most natural experimental mean for detecting WIMPs is based on the mea- 
surement of the recoil that a nucleus of an appropriate detector would suffer, when 
hit by a WIMP. Various methods may be employed for the detection of the recoil 
energy 0) . With the present experimental sensitivities, the only signature for dis- 
entangling the signal from the background is based on the annual-modulation effect 



of the signal, due to the composition of the solar-system speed relative to the dark 
halo with the rotational velocity of the Earth around the Sun 113). As we show 
in the following, the discovery potential of WIMP direct searches confers to this 
experimental mean a prominent role in the investigation of particle dark matter. 

Experiments for WIMP direct detection provide a measurement (or an 
upper bound) of the differential event rate 

-Nr^ fdvmv4^iv,E,), (2) 



dEn mw J dEn 

where Nt is the number of the target nuclei per unit of mass, mw is the WIMP 
mass, pw is the local WIMP matter density, v and f{v) denote the WIMP velocity 
and the velocity distribution function in the Earth frame {v = \v\) and da/dE^ is 
the WIMP-nucleus differential cross section. The nuclear recoil energy, E^, is given 
by En = m^gjf^(l — cos6*)/mN, where 6* is the scattering angle in the WIMP- 
nucleus center-of-mass frame, is the nuclear mass and mred is the WIMP- 
nucleus reduced mass. Eq.(0) refers to the case of a monoatomic detector, like the 
Ge detectors. Its generalization to more general situations, like for instance the case 
of Nal, is straightforward. In what follows, pw will be factorized in terms of the 
local value for the total non-baryonic dark matter density pi and of the fractional 
amount of density, ^, contributed by the candidate WIMP, i.e. pw = i ■ Pi- For pi 
we use the range 0.2 GeV cm~^ < Pi < 0.7 GeV cm~'^, where the upper side of the 
range takes into account the possibility that the matter density distribution is not 
spherical, but is described by an oblate spheroidal distribution HI' HI) . 

The WIMP-nucleus differential cross section may conveniently be split 
into a coherent part and a spin-dependent one = (;^)^ + (^^)^^! whose 



generic features are discussed in the seminal paper of Ref. iD. To compare the- 
oretical expectations with experimental data, and experimental data of different 
detectors among themselves, it is useful to convert the WIMP-nucleus cross-section 
into a WIMP-nucleon cross section. Under the hypothesis that the coherent cross- 
section is dominant and the WIMP couples equally to protons and neutrons (at least 
approximately), the WIMP-nucleus cross section may be expressed in terms of a 
WIMP-nucleon scalar cross section (Jgcaiar""'' ® 



da ^ ( da \ ^ F'^{q) / 1 + mw/mp^ 
dE~R ~ [dE^Jc ~ [l + mw/m^, 



^Vttr\ (3) 



where rrip is the proton mass, A is the nuclear mass number, E"^"-^ is the maximal 
recoil energy and F{q) is the nuclear form factor for coherent interactions, usually 



parametrized in the Helm form \£IV . In the following we assume that the previous 
conditions are satisfied, so that, by using Eq. (H), a WIMP-nucleon scalar cross 
section cricailv"^^ "^^Y be derived from the WIMP-nucleus cross section. 

From the general formula in Eq. (0), it is clear that, in order to extract a 
WIMP-nucleus cross section from the experimental data, one has to use a specific 
expression for the velocity distribution function f{v) (in writing Eq. (0) we have 
already assumed that the WIMP phase-space distribution function is factorizable 



as p(r) ■ f{v), though this is certainly not the most general case The usual 

choice for f{v) is the isotropic Maxwell-Boltzmann distribution in the galactic rest 
frame, as derived from the isothermal-sphere model. The results from WIMP direct 
measurements, as derived by employing the standard isotropic Maxwell-Boltzmann 
distribution, will be presented in the following section; the consequences of using 
other WIMP distributions will be discussed in Sect. 4.1.2. 

4.1.1 Results from WIMP direct measurements 

In Fig. |l] we summarize the most recent experimental results of WIMP direct mea- 



surements 3 cEp e3' l3' e3' B Ezp eH) in terms of as a function of the 

WIMP mass. In deriving the plots of Fig. 0, the isotropic Maxwell-Boltzmann dis- 
tribution has been used for the velocity distribution f{v) and the local WIMP density 
has been set at the default value pi = 0.3 GeV cm~^. The open curves denote upper 
bounds, the closed one denotes the 3-cr region, derived from the annual-modulation 
effect measured by the DAMA Collaboration 



It has been shown S' 0' B HI' H EIP that the DAMA annual-modulation 
effect may be interpreted as due to relic neutralinos, whose relic abundance may also 
be in the range of cosmological interest. Comparisons of the experimental data of 
Ref. HI) with susy calculations have also been presented in Refs. ^ S El' ED . 



In Sect. 5 we will discuss these properties, by proving that the current direct ex- 
periments for WIMPs, when interpreted in terms of relic neutralinos, probe regions 
of the supersymmetric parameter space compatible with all present bounds from 
accelerators. 

4.1.2 Dependence on the WIMP distribution function in the halo 

Let us now discuss some properties related to the use of WIMP distribution functions 
different from the standard Maxwell-Boltzmann distribution. Recent investigations 
have shown that deviations from this standard scheme, either due to a bulk rotation 
of the dark halo E3' El) or to an asymmetry in the WIMP velocity distribution 
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Figure 1: Experimental results of WIMP direct measurements ^s-ss) terms of 
C,(Jscaiai"^^ as a function of the WIMP mass mw- The isotropic Maxwell-Boltzmann 
distribution has been used for the velocity distribution f{v) and the local WIMP 
density has been set at the default value pi = 0.3 GeV cm~^. The open curves 
denote upper bounds, the closed one denotes the 3-a region, derived from the annual- 



modulation effect measured by the DAMA Collaboration 
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influence the determination of tlie WIMP-nucleus cross sections from 
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tfie experimental data in a sizeable way. In Ref. also triaxial matter distribu- 
tions are considered; in the present paper deviation from sphericity in the WIMP 
matter distributions are taken into account only through the physical range allowed 
for the value of pi (see our previous comment on pw after Eq. (||)). In the typical 



plots of cr. 
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mw, as the ones displayed in Fig. |I], the effect introduced by 



the mentioned deviations from the Maxwell-Boltzmann is generically to elongate 
the contours towards larger values of mw- This is for instance the case for the the 
annual-modulation region of the DAMA Collaboration 0) . In Fig. 3 of Ref. 0) it 
is shown that, by implementing the dark halo with a bulk rotation according to the 
treatment in Ref. , the annual-modulation region moves towards larger values of 
the WIMP mass, with an elongation which brings the right-hand extreme from the 
value of ~ 150 GeV to ~ 200 GeV. A similar effect is obtained by introducing an 
asymmetry in the WIMP velocity distribution f{v): Fig. 4 of Ref. 0) illustrates 
this point; notice that this asymmetry effect also pushes somewhat downwards the 



annual-modulation region. We emphasize that all these effects are extremely im- 
portant, when experimental results of WIMP direct detection are being compared 
with theoretical models for specific candidates. This point has been overlooked in 
most analyses in terms of relic neutralinos. 

4.1.3 WIMP local density 

The values to be assigned to the WIMP local density pw = ^-pi have to be consistent 
with the values of the WIMP relic abundance, as derived from evaluations for any 
given specific WIMP candidate. For the case considered in detail in the following, the 
neutralino, we proceed in the following way. The relic abundance Vtyji^ is evaluated 
in specific supersymmetric schemes (see Sect. V). When fi^/i^ > (i7m^^)mm, where 
{^mh'^)min is the minimum value of VLmh? compatible with halo properties, we simply 
set p^ = pi (i.e., ^ = 1). When fl^h"^ < {Qmh'^)min, the neutralino cannot be 
the unique cold dark matter particle, thus we assign to the neutralino a resettled 
local density p^ = pi x Vt^h^ / {yLmh'^)min (i-e., ^ = Q^h"^ / {Qrnh'^)min) Thus, 
summarizing, 



Px = Ph when Q^h^ > {Qrnh'^)min (4) 

Px = ^^^^ ^X^^ < (^m/i^)mm- (5) 

\\LifYill )min 

We stress that in our analyses we consider various scenarios, where neu- 
tralinos may constitute the dominant cold dark matter component or only a small 
fraction of it. As discussed below (see Sect. 4.3), in view of the properties related 
to detectability of relic neutralinos, it is very important not to disregard neutralino 
configurations with small relic abundance, where rescaling applies. 

4.1.4 Sensitivity range for current WIMP direct searches 

In the present paper we focus our analysis to the WIMP mass range which, in the 
light of the experimental data summarized above and of the previous considerations 
on the astrophysical uncertainties, appears particularly appealing: 

40 GeV <mw< 200 GeV. (6) 

Notice that the mass range of Eq. (^ is quite appropriate for neutralinos. 
Actually, the lower extreme is indicative of the LEP lower bound on the neutralino 
mass (in the calculations performed in the present work the actual lower bound 



for m^, dependent on the other susy parameters, is employed, according to the 
constraints given in 0)). As for the upper extreme, we notice that, though a generic 
range for might extend up to about 1 TeV, requirements of no excessive fine- 
tuning would actually favour an upper bound of order 200 GeV, in accordance 
with Eq. (§). 

In what follows we will discuss the discovery potential of WIMP direct 
searches for WIMPs in the mass range of Eq.(^. Particular attention will be paid 
to capabilities of the present experiments; according to the discussion previously 
made, their sensitivity range, in case of WIMPs whose coherent interactions with 
ordinary matter are dominant over the the spin-dependent ones, may be stated, in 
terms of the quantity '^scalar 

4 ■ 10-1° nbarn < ^a^r^ < 2 ■ 10"^ nbarn. (7) 

We will hereafter refer to region R as the one in the space mw ~ ^cricailr"^^ 
which is defined by Eqs. The region R represents the sensitivity region already 

under exploration with present detectors. 

Our analysis, based on an interpretation of experimental data in terms 
of relic neutralinos, will show by how much the WIMP direct searches probe the 
supersymmetric parameter space. 

4.2 WIMP indirect searches 

Indirect searches for WIMPs aim at the measurements of signals which fall essentially 
into two categories (for a review and relevant references see, for instance, Ref. 0)): 

i) Signals due to WIMP pair-annihilations taking place in our galactic 
halo. The most interesting outputs of these annihilation processes are: a) fluxes of 
neutrinos and gammas (which might be disentangled from the background in case of 
halo dumpiness), b) a narrow gamma-gamma line, c) exotic components in cosmic 
rays: antiprotons, positrons, antideuterons. 

ii) Signals due to WIMP pair-annihilations taking place in the interior of 
celestial bodies. These signals would consist of fluxes of up-going muons in a neu- 
trino telescope, generated by neutrinos which are produced by the pair annihilations 
of WIMPs captured and accumulated inside the Earth and the Sun. 

For the discussion to follow, it is important to notice that in case (i) the 
detection rates are proportional to ■ aann, whereas in case (ii) the detection rates 
are proportional to pw ■ cricailr"^^ (^^ iact, in this latter case, the annihilation rate 



inside the celestial body is proportional to the capture rate of the WIMPs by the 
body, which is in turn proportional to pw ■ (J^^^yw"^"^)- 



4.3 General arguments on detectability of WIMPS versus their local density and 
average relic abundance 

One particular type of relic particle (say, a WIMP) may constitute a dominant or a 
subdominant dark matter candidate. Also the latter case would be of great physical 
interest, provided this relic particle has some chance to be detected. Thus, a crucial 
question is: How do experimental detectablities of WIMPs depend on the WIMP 
(local and cosmological) densities ? 

This question may be answered by simple general arguments which we 
already presented in a number of papers ^ " EU' and which we 

briefly review here. Our arguments are exemplified by properties of relic neutralinos, 
in case these particles decouple from the primordial plasma when in thermodynamic 
equilibrium. It is obvious that these same arguments apply equally well to other 
realizations of WIMPs, under the same decoupling conditions. 

We consider the rates of direct signals and of signals due to pair annihilation 
in celestial bodies, separately from the rates of signals due to neutralino-neutralino 
annihilations taking place in the galactic halo, because of their different dependence 
on local density and cross sections. 

4.3.1 Rates of direct detection and of signals due to pair annihilation in celestial 
bodies 

As was seen above, in both of these experimental measurements the detection rate R 
is proportional to the product ■ (J^cailr"^^ ■ Thus, taking into account the rescahng 
properties of p^ (Eqs. (H)-(||)), we have that R behaves as follows 



R ^ ^^c^lar°''^ when n^h"^ > {Vtmh^)^in (8) 

O 7,2 (nucleoli) 

R oc ''f attr^ oc ^.s^, when U^h' < {Q^h'Un. (9) 

Now, the cross sections and o"ann, as functions of any generic 

coupling parameter behave similarly {i.e. they usually both decrease or increase 
in terms of variations of that parameter), because of crossing symmetry. Thus, 
for instance, in supersymmetric schemes where the neutralino is the LPS, cTg^i^^^''""'* 
and (Jann are both increasing functions of tan (3, when the relevant processes are 
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Figure 2: Qualitative behaviour of the rates of direct detection and of signals due 
to pair annihilation in celestial bodies as functions of a generic coupling parameter 
( (C,r denotes the value of ( at which ^^^^ = {^mh'^)min)- Notice the correlation 
between the rates R and the neutralino relic abundance Vt^h'^ . 



mediated by Higgs bosons. Usually, CTgcaiar'"'^'' increases somewhat faster than Oann, or 
approximately at the same rate. Since, at the same time, Vt^h^ oc (< aannV >)^^, the 
typical behaviour of the rate R (for the processes under study) is the one displayed 
in Fig. 1^. For small values of C, both cTgcaiar""^ ^^'^ '^ann are small (and then Vt-Ji^ 
is large) and R grows proportionally to o'gcaiar""^'' • the strength of the coupling 
increases, in the region beyond the value Cr (at which Vty^h"^ = (fim^^)mm), Eq. 
(0) applies: the rate R still increases (though less rapidly than before rescaling), 
or remains approximately flat. Therefore, as far as direct detection and indirect 
detection through pair annihilation in celestial bodies are concerned, we obtain that 
the detectability of relic neutralinos is usually favoured for neutralinos of small Q^h"^ , 
that is for neutralinos which comprise only a subdominant dark matter component. 

We stress that, due to this last property, in the numerical analyses of the 
properties of relic neutralinos, it is important to keep into consideration also susy 
configurations where rescaling applies. In fact, by neglecting these configurations, 
one would disregard relic neutralinos which have the largest chances of being de- 
tected. Because of this fact, at variance with most analyses by other authors, our 
scanning of the susy parameter space always included configurations entailing small 
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Figure 3: Qualitative behaviour of the rates of signals due to neutralino-neutralino 
annihilations taking place in the galactic halo as functions of a generic coupling 
parameter ( ((r denotes the value of ( at which fl^h'^ = {flrnh'^)min)- Notice the 
correlation between the rates R and the neutralino relic abundance fl^h"^. 



neutralino relic abundances (see, for instance, 

The numerical evaluations, which we present in the next section, show 
that the actual situation is even better; in fact, it turns out that a number of 
neutralino configurations of cosmological interest, though disfavoured by the previous 
arguments, may reach the level of detectability by current experiments. Actually, we 
will prove that the annual-modulation effect, measured by the DAMA Collaboration, 
is compatible with an interpretation in terms of neutralinos, part of which with relic 



abundance of cosmological interest. Due to the connection between a. 
Qyh'^, this property is far from being trivial. 
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and 



4.3.2 Rates of signals due to neutralino-neutralino annihilations taking place in 
the galactic halo. 



At variance with the previous case, the detection rates have now a quadratic de- 
pendence on the neutralino density in the halo, R (x ■ aann', thus, by using Eqs. 
©"(i)' ^6 have the following behaviour for R 



R when ily^h > {Vtmh )mm (10) 

R OC ijTr-J^ y CTann OC when n^h'^ < {^lmh'^)min- (H) 

Therefore, in this the generic couphng parameter C, increases, the 

detection rate first increases, as long as the couplings are small; afterwards, for 
values of C larger than C^, when rescaling applies, the detection rate decreases (see 
Fig. I). 

It follows that for processes depending on pair-annihilation in the halo the 
maximal rates occur for values of the relic abundance around the value {VLmh'^)min- 
Subdominant neutralinos are disfavoured for detectability by this type of signals as 
compared to neutralinos with a relic abundance around the value ^^^'^ — {^mh^)min- 

Finally, we notice that all the arguments reviewed in this section were 
already presented in our previous papers (see, for instance, m 0)) , further 
corroborated by extensive numerical evaluations S S' HO). Amazingly, in a 



recent preprint ED , where considerations similar to ours are repeated, all previous 
references have been overlooked. 



5 Supersymmetric dark matter 

We turn now to calculations performed in specific supersymmetric schemes, assum- 
ing that R-parity is conserved, and thus that the LSP is stable. The nature of the 
LSP depends on the susy-breaking mechanism and on the specific regions of the 
susy parameter space. We consider here gravity-mediated schemes, and domains 
of the parameter space, where the LSP is the neutralino. Extensive calculations on 
relic neutralino phenomenology in gravity-mediated models have been performed 
(see, for instance, Refs.3i-4i,49,5i-55,57,59-65))^ 

Our analysis will follow the one reported in Ref. EO) and will show by how 
much the WIMP direct searches probe the supersymmetric parameter space. We 
remark that, in the case of neutralinos, the assumption about the dominance of the 
coherent cross section over the spin-dependent one is, in general, largely satisfied, 
except for values of o'f^aiar'"'^'' which are far below the present experimental reach 

HD. 



5.1 Susy models 

The analysis of Ref. 0) was performed in the Minimal Supersymmetric extension of 
the Standard Model (MSSM) in a variety of different schemes, from those based on 



universal or non-universal supergravity, with susy parameters defined at the grand 
unification scale, to an effective supersymmetric model defined at the Electro- Weak 
(EW) scale. Here we only report about results in universal supergravity and in the 
effective scheme at EW scale; for the other cases and further details we refer to 
Ref. 0). 

The essential elements of the MSSM are described by a Yang-Mills La- 
grangian, the superpotential, which contains all the Yukawa interactions between 
the standard and supersymmetric fields, and by the soft-breaking Lagrangian, which 
models the breaking of supersymmetry. The Yukawa interactions are described by 
the parameters h, which are related to the masses of the standard fermions by the 
usual expressions, e.g., rrit = 2, riif, = h^vi, where Vi are the few's of the two Higgs 
fields. Hi and H2- Implementation of this model within a supergravity scheme leads 
naturally to a set of unification assumptions at a Grand Unification (GUT) scale, 
Mqut'- 

i) Unification of the gaugino masses: Mi^Mcur) = fni/2, 

ii) Universality of the scalar masses with a common mass denoted by mo: 
mi{MGUT) 

= mo, 

iii) Universality of the trilinear scalar couplings: A\Mgut) = A''-{Mgut) = 

A"(MGt/T) 

This scheme will be denoted here as universal SUGRA (or simply SUGRA). 
The relevant parameters of the model at the electro-weak (EW) scale are obtained 
from their corresponding values at the Mqut scale by running these down according 
to the renormalization group equations (RGE). By requiring that the electroweak 
symmetry breaking is induced radiatively by the soft supersymmetry breaking, one 
finally reduces the model parameters to five: mi/2, mo, Aq, tan/5(= V2/V1) and sign 
/i. Here the parameters are varied in the following ranges: 50 GeV < mi/2 < 
1 TeV, mo < 1 TeV, —3 < Aq < +3, 1 < tan (3 < 50. Notice that a common upper 
extreme for the mass parameters has been used, and generically set at the value of 1 
TeV, as a typical scale beyond which the main attractive features of supersymmetry 
fade away. However, fine-tuning arguments actually set different bounds for mo and 
mi/2 (in universal SUGRA and in non-universal SUGRA) EIP: mi/2 < hundreds 
of GeV, whereas mo < 2 — 3 TeV. The mg ~ 2 — 3 TeV window, which is not 
specifically included in the results reported here, will be analysed in a forthcoming 

properties of relic neutralinos in this large mo regime have 
been analyzed. We remark that the phenomenology of relic neutralinos is also very 
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sensitive to other parameters, such as the top quark mass rrit and the strong couphng 
as- For these parameters, we use here their 95% CL ranges: nit = (175 ± 10) GeV 
and as{Mz) = 0.118 ± 0.004. 

Models with unification conditions at the GUT scale represent an appealing 
scenario; however, some of the assumptions listed above, particularly ii) and iii), are 
not very solid, since, as was already emphasized some time ago universality 
might occur at a scale higher than Mqut ~ 10^^ GeV, e.g., at the Planck scale. 
More recently, the possibility that the initial scale for the RGE running, M/, might 
be smaller than Mqut ~ 10^^ has been raised HD, on the basis of a number 
of string models (see for instance the references quoted in HD). In Ref. 0) it is 
stressed that M/ might be anywhere between the EW scale and the Planck scale, 
with significant consequences for the size of the neutralino-nucleon cross section. 

An empirical way of taking into account the uncertainty in Mj consists in 
allowing deviations in the unification conditions at Mqut- For instance, deviations 
from universality in the scalar masses at Mqut, which split from may be 
parametrized as 

MKMgut) = mlil + 6,) (12) 

This is the case of non-universal SUGRA (nuSUGRA) that we considered 
Refs. 23' 0^ . Further extensions of deviations from universality in SUGRA 



m 



models which include squark and/or gaugino masses are discussed, for instance, in 
^ 0) In this note we do not report results in nuSUGRA schemes; we refer to 



Ref. iO) for discussions on this case. 

The large uncertainties involved in the choice of the scale Mj make the 
SUGRA schemes somewhat problematic: the originally appealing feature of a uni- 
versal SUGRA with few parameters fails, because of the need to take into considera- 
tion the variability of Mj or, alternatively, to add new parameters which quantify the 
various deviation effects from universality at the GUT scale. It appears more con- 
venient to work with a phenomenological susy model whose parameters are defined 
directly at the electroweak scale. We denote here this effective scheme of MSSM by 
effMSSM. This provides, at the EW scale, a model, defined in terms of a minimum 
number of parameters: only those necessary to shape the essentials of the theoret- 
ical structure of an MSSM, and of its particle content. Once all experimental and 
theoretical constraints are implemented in this effMSSM model, one may investigate 
its compatibility with SUGRA schemes at the desired Mj. 

In the effMSSM scheme we consider here, we impose a set of assumptions 



at the electroweak scale: a) all trilinear parameters are set to zero except those of 
the third family, which are unified to a common value A; b) all squark soft-mass 
parameters are taken degenerate: rriq- = 171^; c) all slepton soft-mass parameters are 
taken degenerate: m^- = mf d) the U{1) and SU{2) gaugino masses, Mi and M2, 
are assumed to be linked by the usual relation Mi = (5/3) tan^ 6*14/ M2 (this is the 
only GUT-induced relation we are using, since gaugino mass unification appears 
to be better motivated than scalar masses universality). As a consequence, the 
supersymmetric parameter space consists of seven independent parameters. We 
choose them to be: M2, /i, tan/?, m^, m^, m^-, A and vary these parameters in the 
following ranges: 50 GeV < M2 < 1 TeV, 50 GeV < < I TeV, 80 GeV < tua < 
1 TeV, 100 GeV < mq,mi < 1 TeV, -3 < A < +3, 1 < tan/3 < 50 (m^ is the 
mass of the CP-odd neutral Higgs boson). 

The effMSSM scheme proves very manageable for the susy phenomenology 
at the EW scale; as such, it has been frequently used in the literature in connection 
with relic neutralinos (often with the further assumption of slepton/squark mass 
degeneracy: rriq = mj) S S S S §3' 0). Notice that we are not assuming 
here slepton/squark mass degeneracy. 

We recall that even much larger extensions of the supersymmetric models 
could be envisaged: for instance, non-unification of the gaugino masses 1^' 
and schemes with CP-violating phases . Here we limit our considerations to the 
two schemes previously defined: universal SUGRA and effMSSM; for nuSUGRA we 
refer to Ref. © . 

The neutralino is defined as the lowest-mass linear superposition of photino 
(7), zino {Z) and the two higgsino states (ifi, iJ|): x = (^il + (^2Z + a^H^ + 04-^2- 
Hereafter, the nature of the neutralino is classified in terms of a parameter P, defined 
as P = af + a^. The neutralino is called a gaugino when P > 0.9, a higgsino when 
P < 0.1, mixed otherwise. 

For more details concerning theoretical aspects involved in our calculations 
we refer to Refs. 23' © 0). Accelerators data on supersymmetric and Higgs bo- 
son searches (CERN e~^e~ collider LEP2 and Collider Detector CDF at Fermilab) 
provide now rather stringent bounds on supersymmetric parameters. CDF bounds 
are taken from The new LEP2 bounds are taken from Erp \L2i' . 

As compared to the calculations presented in Ref. 0), we have now im- 
plemented effects due to those radiative corrections to the couplings of the neutral 
Higgs bosons to the quarks which may be sizeable at large tan (3 . These radia- 
tive corrections affect the calculation of the neutralino-nucleus cross section and of 
the neutralino cosmological relic abundance. Notice that the correction to the rela- 



tion between the b quark mass and its Yukawa couphng, imphed by these radiative 
corrections, enters also in the calculations of the 6 — s + 7 decay For the 

SUGRA model discussed above, it affects also the boundary conditions at the GUT 
scale for the b Yukawa coupling EH). This in turn affects the radiative symmetry 



breaking mechanism and the low-energy Higgs and sfermion spectra . All these 
effects are included in our calculations. 

Finally, we notice that a new experimental constraint on supersymmetric 
parameters may also be derived from the recent accurate experimental determination 
of the muon anomalous magnetic moment H) ; this measurement provides the value 
a^(exp) = 11659202(14)(6) x 10-^°(1.3 ppm), where = {g - 2)/2. This data, 
if compared with the theoretical evaluations in Ref. H), would show a deviation 
of 2.6 a from the standard-model prediction. This has determined an outburst of 
theoretical papers HHP , where this possible deviation is attributed to supersymmetry, 
and the relevant implications derived. However, other more recent standard-model 
evaluations of a,, are in better agreement with the experimental data of 

Ref. u) (see also Refs. 0, 0) for detailed discussions of the standard-model 



calculations of a^). Thus, for the time being, it appears safer to use the data of Ref. 

as a constraint on susy, rather than a sign of it. Employing the theoretical results 
of Refs. 1^' 0' J the contribution of supersymmetry to the anomalous moment 
is constrained by —140 < a^"*^ ■ 10^^ < 890 . This constraint has been implemented 
in our present scanning of the supersymmetric parameter space. 

5.2 Numerical results 

We turn now to the presentation of our results. In Figs. we give the scatter 
plots for (jI^^yIt"^^ versus ^^^'^ different schemes: universal SUGRA, and 

effMSSM. For the SUGRA scheme we only display the results corresponding to 
positive values of fi, since, for negative values, the constraint on 6 — > s + 7 implies 
a large suppression of c^l^^i^^"^^ ■ The calculations of have been performed 

with the formulae reported in Refs. 0); set 1 for the quantities rriq < qq >'s has 
been used (see Ref. 0) for definitions); the evaluation of Vt^K^ follows the procedure 
given in 0) . The two horizontal lines bracket the sensitivity region defined by Eq. 
(^, when ^ is set equal to one. The two vertical lines denote the favorite range for 
VLmh'^, 0.05 < VLmh"^ < 0.3. 

Figs. provide a first relevant result of our analysis: the present experi- 
mental sensitivity in WIMP direct searches allows the exploration of supersymmetric 
configurations compatible with current accelerator bounds. A number of configura- 
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Figure 4: Scatter plot of al^!^^lr°^^ versus Q^h'^ for universal SUGRA. Set 1 for the 

quantities ruq < qq > 's is employed 1^ . Only configurations with positive /i are 
shown and is taken in the range of Eq. (j^. The two horizontal lines bracket the 
sensitivity region defined by Eq. (Qj, for ^ = 1. The two vertical lines denote the 
range 0.05 < flrnh"^ < 0.3. The region above Q^h"^ = 0.7 is excluded by current limits 
on the age of the Universe. Dots (crosses) denote gaugino (mixed) configurations. 



tions stay inside the region of cosmological interest, also in the constrained SUGRA 
scheme. The region of experimental sensitivity and cosmological interest is covered 
with an increasingly larger variety of supersymmetric configurations as one moves 
from SUGRA to effMSSM; this latter fact is expected from the intrinsic features of 
the various schemes. 

Once a measurement of the quantity ■ (t^^^w'^^'' performed, values for 
the local density versus the relic abundance Vt-^h'^ may be deduced by proceeding 
in the following way 0) : 

-, \ . 1 i 1 r (nucleon) 1 / (nucleon) i r (nucleon) i i , 

1) IS evaluated as [p^ ■ '\expt / t^scaiar ^ where [p^ ■ al^^^.^^ '\expt denotes 

the experimental value, and CTgcaiar""^^ calculated as indicated above; 2) to each 
value of p^ one associates the corresponding calculated value of fl^h"^. The scatter 
plot in Fig. ^ is derived from the lowest value of the annual-modulation region 
of Ref. [p^/(0.3 GeV cm~'^) ■cTgcaiar'"'^'' ]expi = 1 ■ 10~^ nbarn, and by taking 
in the range of Eq. (B). This plot, obtained in case of effMSSM, shows that the 
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X 

Figure 5: Scatter plot of a^^!^^^^°^^ versus ^^^^ for effMSSM. Notations as in Fig. 0. 
Dots denote gauginos, circles denote higgsinos and crosses denote mixed configura- 
tions. Both signs of fi are shown. 



most interesting region, i.e. the one with 0.2 GeV cm < — O-"^ GeV cm ^ 
and 0.05 < Vlmh^ < 0.3 (cross-hatched region in the figure), is covered by susy 
configurations probed by the WIMP direct detection. 

Let us examine the various sectors of Fig. |^. Configurations above the 
upper horizontal fine are incompatible with the upper limit on the local density 
of dark matter in our Galaxy and must be disregarded. Configurations above the 
upper slanted dot-dashed line and below the upper horizontal solid line would im- 
ply a stronger clustering of neutralinos in our halo as compared to their average 
distribution in the Universe. This situation may be considered unlikely, since in 
this case neutralinos could fulfill the experimental range for p^, but they would 
contribute only a small fraction to the cosmological cold dark matter content. For 
configurations which fall inside the band delimited by the slanted dot-dashed lines 
and simply-hatched in the figure, the neutralino would provide only a fraction of 
the cold dark matter at the level of local density and of the average relic abun- 
dance, a situation which would be possible, for instance, if the neutralino is not 
the unique cold dark matter particle component. To neutralinos belonging to these 
configurations one should assign a resca/ed local density (see Sect. 4.1.3). 



. Bottino, F, Doneto, N, Fornonoo, S, Scopol (2001) 



10-8 
10-- 























'<:* • 




























X 


■.•••■■1 


'/ ' 1 






^ . r ; ■ 








J. 

■ .f 






6 X , , X 


; '1 






Qx ^ X 


, ,1,,,,, 







Q h2 



Figure 6: Scatter plot of versus ^^^"^ /^'^ ejJMSSM. This plot is derived from 



the experimental value \p^/(0.3 GeV cm ^) -cr. 



(nuclcon)-i 
scalar 



J expt 



1 ■ 10 ^ nharn and by 



taking in the range of Eq. according to the procedure outlined in the text. Set 
1 for the quantities rUq < qq > 's is employed ^E) . The two horizontal lines delimit 
the range 0.2 GeV cm~^ < Px — ^■'^ GeV cm~^; the two vertical ones delimit the 
range 0.05 < ^mh"^ ^ 0.3. The region above ^^^^ ~ ^-^ ^■^ excluded by current limits 
on the age of the Universe. The band delimited by the two slanted dot-dashed lines 
and simply hatched is the region where rescaling of pi applies. Dots denote gauginos, 
circles denote higgsinos and crosses denote mixed configurations. 



We remind that the scatter plot in Fig. ^ refers to a representative value 
of [p^ ■o'icaiar'"'^''] ^i^sidc the currcnt experimental sensitivity region, thus the plot in 
Fig. P shows that current experiments of WIMP direct detection are probing relic 
neutralinos which may reach values of cosmological interest, but also neutralinos 
whose local and cosmological densities may provide only a very small fraction of these 
densities. These properties were anticipated by the general arguments previously 
discussed in Sect. IV. C. 

For sake of comparison with specific experimental results, we provide in 
Figs. the scatter plots for the quantity ^cic^iar""^ versus in the two su- 
persymmetric schemes. In these figures the solid line denotes the frontier of the 
3a annual-modulation region of Ref. 0) , when only the uncertainties in pi and in 
the dispersion velocity of a Maxwell-Boltzmann distribution, but not the ones in 
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Figure 7: Scatter plot of ^cr^cailr"^^ versus m^^ in case of universal SUGRA. Set 
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1 for the quantities niq < qq >'s is employed Eli), Crosses (dots) denote con- 
figurations with f^x^^ ^ 0-0^ (VlJt^ < 0.05j. The solid contour denotes the 3a 



annual-modulation region of Ref. 1^ (with the specifications given in the text). 



other astrophysical quantities, are taken into account. As discussed in the Intro- 
duction, effects due to a possible built rotation of the dark halo or to an asymmetry 
in the WIMP velocity distribution would move this boundary towards higher val- 



ues of rriy^. Our results in Figs, 
annual-modulation region of Ref. 



show that the susy scatter plots reach up the 
' , even with the current stringent bounds from 
accelerators (obviously, more easily in effMSSM than in SUGRA). 

Finally, we recall that use of set 2 for the quantities ruq < qq >'s instead 
of set 1 (see, for definitions, Ref. ) would entail an increase of about a factor 3 
in all the scatter plots of Figs. 



6 Conclusions 



In this work, after a short general introduction, we have shown that the current 
direct experiments for WIMPs, when interpreted in terms of relic neutralinos, are 
indeed probing regions of the supersymmetric parameter space compatible with 
all present bounds from accelerators. We have quantified the extent of the explo- 
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Figure 8: Same as in Fig. in case of ejfMSSM. 



ration attainable by WIMP direct experiments in terms of different supersymmetric 
schemes, from a SUGRA scheme with unification assumptions at the grand unifi- 
cation scale to an effective model, effMSSM, at the electroweak scale. It has been 
stressed that, due the large uncertainties in the unification assumptions in SUGRA 
schemes, the effMSSM framework turns out to be the most convenient model for 
neutralino phenomenology. 

We have proved that part of the configurations probed by current WIMP 
experiments entail relic neutralinos of cosmological interest, and, a fortiori, also neu- 
tralinos which might contribute only partially to the required amount of dark matter 
in the Universe. This last property was anticipated by the arguments presented in 
Sect. 4.3 and confirmed by our numerical results. The cosmological properties have 
been displayed in terms of a plot of the local density versus the average relic abun- 
dance, i.e. in a representation which proves particularly useful to summarize the 
properties of relic neutrahnos (see Fig. ||). 

The question: Are relic neutralinos of very low (local and cosmological) 
densities detectable by current experiments of WIMP direct detection ? finds a 
straightforward and affirmative answer in the versus plot. Direct detectabil- 
ity is possible even for neutralino densities quite minuscule as compared to the ones 
of cosmological interest. 



In the present note, our discussions were mainly focussed on implications 
of direct detection results. However, as mentioned in Sect. 4.2, also WIMP in- 
direct searches are quite important. Measurements of up-going muon fluxes from 
the center of the Earth and from the Sun can potentially either find a signal or, 
at least, place significant constraints (though some uncertainties about a possible 
solar-bound population have still to be resolved) Also measurements of low- 

energy antiprotons in space may provide interesting constraints on the susy model 
parameters 1^), in view of the new refined calculations of secondary antiprotons 
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We notice that, by the arguments presented in Sect. 4.3, similarly to the 
case of direct detection, the detectability of relic neutralinos by measurements of 
up-going muon fluxes from the center of the Earth and from the Sun at neutrino 
telescopes is more favourable in case of low local and average densities. 

We have shown that the annual-modulation effect measured by the DAMA 
Collaboration may be interpreted as due to relic neutralinos, which are compati- 
ble with all current constraints from accelerator measurements and WIMP indirect 
searches. 

In our evaluations we have taken into account that the determination of 
the actual sensitivity region in terms of the WIMP-nucleon cross section and of the 
WIMP mass from the experimental data depends quite sizeably on uncertainties of 
various origins, mainly: i) possible effects due to a halo bulk rotation and/or to 
asymmetries in the WIMP velocities distribution, and ii) significant uncertainties 
in the determination of Higgs-quark-quark and neutralino-quark-squark couplings. 
We stress that all these effects have to be taken properly into account, when con- 
clusions about comparison of theory with experiments are drawn. 

Finally, we wish to point out that a susy Higgs boson at a mass of about 
115 GeV, as possibly hinted by the Higgs LEP experiments would fit remarkably 
well in the scenario depicted above . 
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